Stochastic perturbations in open chaotic systems: random versus noisy maps 
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We investigate the effects of random perturbations on fully chaotic open systems. Perturbations can be ap- 
plied to each trajectory independently (white noise) or simultaneously to all trajectories (random map). We com- 
pare these two scenarios by generalizing the theory of open chaotic systems and introducing a time-dependent 
conditionally-map-invariant measure. For the same perturbation strength we show that the escape rate of the 
random map is always larger than that of the noisy map. In random maps we show that the escape rate k and 
dimensions D of the relevant fractal sets often depend nonmonotonic ally on the intensity of the random pertur- 
bation. We discuss the accuracy (bias) and precision (variance) of finite-size estimators of k and D, and show 
that the improvement of the precision of the estimations with the number of trajectories N is extremely slow 
(oc 1/ In N). We also argue that the finite-size D estimators are typically biased. General theoretical results are 
combined with analytical calculations and numerical simulations in area-preserving baker maps. 



PACS numbers: 05.45.-a,05.40.-a 

I. INTRODUCTION 

External perturbations affect almost any observation to be 
made and are usually modeled through simple stochastic pro- 
cesses [ 1 1. In this paper we are interested in stochastic per- 
turbations in open chaotic maps, i.e., discrete-time systems 
which exhibit a transiently chaotic dynamics. Such systems 
appear in a variety of physical situations (scattering, planetary 
astronomy, chemical reactions, fluid dynamics, environmen- 
tal sciences, etc.). In an ensemble-based framework, there are 
two different ways of introducing perturbations, for which the 
following common terminology applies [2|: 

Noisy map: perturbations are applied independently to 
each trajectory. 

Random map: the same perturbation is applied to all 
trajectories simultaneously. 

In fluid advection noisy maps would model the effect of 
molecular diffusion, and random maps would model a ran- 
domly changing velocity field affecting all tracers simultane- 
ously [3-8]. Another physical situation where the random 
map picture applies is wave front propagation through ran- 
domly structured media, e.g. in underwater acoustics l9l [l0l . 

The above distinction can be motivated also from a pre- 
dictability point of view, whereby noisy maps describe models 
with uncertainties, and random maps describe models which 
are perturbed by an a priori known process. Here we are con- 
cerned with the predictability of a typical trajectory, which 
is arbitrarily chosen from an ensemble, and the measures of 
predictability will be defined as averages over this ensemble. 
In the well-studied case of dissipative closed systems the en- 
semble at any time t is taken to be constituted by trajectories 
which are arbitrarily initialized in the infinitely distant past, 
t — » — oo. In the random map framework this ensemble is 
referred to as a random or snapshot attractor, which is a frac- 
tal set if the trajectories are chaotic [2|. A remarkable prop- 
erty of the snapshot attractor is that its geometry and the mea- 
sure supported by it are changing continuously in time, but 



its fractal dimension is constant ||2][6][TT). However, e.g. the 
finite-time average maximal Lyapunov exponent, quantifying 
the finite-lead-time predictability of the typical trajectory, is 
time-dependent [12|. On the other hand, in the noisy map 
framework the average point-wise prediction error at some 
finite lead time cannot be arbitrarily reduced by improving 
the precision of the initial conditions. Thus the random map 
picture may seem to be in stark contrast to the noisy dynam- 
ics, however, it has been shown that fractal snapshot attrac- 
tors constitute building blocks of the noisy stationary attrac- 
tor na. 

In autonomous chaotic open systems the density of trajec- 
tories surviving for a long time inside the system, not leav- 
ing a window of observation, decays exponentially, and it is 
distributed according to the so-called conditionally invariant 
measure |[T3T[T5l . In these systems, predictability concerns 
whether we can foretell the route of escape of a trajectory 
when there are a number of options for that. A measure of 
this predictability of the typical, i.e., arbitrarily initialized tra- 
jectory in the observational window is the uncertainty expo- 
nent. The latter is the scaling exponent of the fraction of pre- 
dictable trajectories with respect to the precision of the ini- 
tial conditions fl6l . For random maps the uncertainty expo- 
nent is thus trivially related to the fractal dimension of the 
boundary of basins from which the trajectories escape through 
different routes. The fractal scaling of the basin boundary is 
time-dependent at finite, practically accessible scales, and so 
predictability (or the 'rate' at which it can be improved) is 
also time-dependent. For noisy maps the basin boundary is 
space-filling, and so, similarly as in case of closed systems, 
predictability cannot be improved at all. 

In this paper we concentrate on simple (fully hyperbolic) 
chaotic open systems and simple stochastic perturbation pro- 
cesses, and focus on the effects of the perturbation strength 
and on the comparison between noisy and random maps. We 
combine the concepts of conditionally invariant [ 15 1 and time- 
dependent but map-invariant ifTSl [TTl [l8l measures, and ar- 
gue that the trajectories in the random open maps are dis- 
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tributed according to a time-dependent conditionally-map- 
invariant measure (to be clarified below). Based on this for- 
malism we obtain that, for the same stochastic perturbation 
process, the escape rate n of the random map is always smaller 
than that of the noisy map. We also investigate the dimen- 
sion D of the relevant fractal sets of random open maps, and 
we discuss the accuracy and precision of finite-size estimators 
of k and D. Under conditions when noise-enhanced trapping 
is observed, i.e., when the average life-time of trajectories is 
'constructively' increased by noise Ifl9ll20l . we find numeri- 
cal evidence that both k and D of the associated random map 
also show a nonmonotonic dependence on the perturbation 
strength. 

The paper is organized as follows. In Sec. [TT] the theory 
of open systems is generalized for random maps, introducing 
a time-dependent version of the conditionally invariant mea- 
sure. In Sec. Ill we provide general relations for the escape 



rate in autonomous, random, and noisy maps, followed by the 
discussion of finite-size estimators. The analogous investiga- 
tion for the fractal dimension appears in Sec. |IV] Finally, our 
main conclusions are summarized in Sec. [V] 



II. MEASURES OF OPEN MAPS 

We consider the temporal evolution over t = 0, . . . , T of 
an ensemble of n = 1, • • • , N trajectories under the action 
of the map Xt+i — ft(xt) in a d-dimensional phase space 
x E X, when on each iteration with respect to t the map f t is 
chosen from an ensemble according to some probability dis- 
tribution (6). We assume that members of this ensemble are 
invertible, i.e., xt — ff (xt+i), open, and fully chaotic (to 
be clarified below). Equivalently, we can say that the mapping 
rule / depends on a control parameter a, and we consider per- 
turbations around a fixed value a* as a = a t = a* +<5£t, where 
5 is the strength of the perturbation, and £ t 's are independent 
identically distributed random variables (e.g. Gaussian with 
zero mean and unit variance), which in general vary across 
different trajectories n (but do not depend on x n ). Altogether, 
the dynamics of the ensemble is written as: 



x t+i, 



ft,n(x t ,n) = /(#t,nj 



0) 



According the theory of open maps |[T5ll22l . for t — > ±oo al- 
most every trajectory leave a finite region of the phase space 
r C X in which they exhibit some nontrivial dynamics. A 
central quantity in our analysis will be the probability density 
function p(x, t) of surviving trajectories in x E r up to time t, 
which is obtained by dividing the number of trajectories in an 
e-neighborhood of x by the total number of surviving trajec- 
tories N(t) (in the limit of e -> when N(0) -> oo). While 
N(t)/N(0) — > for t — > oo, the normalized density p(x,t) 
may approach a nontrivial density and be used to define a mea- 
sure (dp = pdVx, where dVx denotes a phase space volume 
element in X). Next we discuss in detail the properties of 
this measure in the cases of autonomous-, random-, and noisy 
maps. 



A. Autonomous maps 

In the unperturbed case, 5 = in Eq. ([!]), the map is au- 
tonomous, and the following results are known from tran- 
sient chaos theory fl5l : the dynamics is governed by a time- 
invariant nonattracting chaotic set in r, also called a chaotic 
saddle, which is composed of the points that do not leave 
r under the action of the mapping ([T]i in either direction 
t — > ±oo (23). For fully chaotic maps this is a zero measure 
fractal set composed by the intersection of the set of points 
that never leave r for t — > oo (stable manifold) with the set 
of points that never leave r for t —> — oo (unstable mani- 
fold). The normalized density p(x, t) converges to a well de- 
fined stationary density p(x) for t — > oo. The measure p 
associated with p(x) is said to be conditionally invariant (in 
brief c-measure and c-density, respectively), because for any 
set A C r it obeys the following relation |[T3l[T4l : 



p(A) 



(2) 



where /~ 1 (_T) C r corresponds to the set of points that 
do not escape r over one iteration of /. Because of a 
constant rescaling given by the denominator in Eq. (|2]), c- 
measures are not invariant under the map /, i.e., not /- 
or map-invariant (TH). For clarity, we can refer to them 
as conditionally-map-invariant. C-measures of autonomous 
maps are time-invariant, however. The c-measure associated 
with p(x) is a probability measure [18], indicating the chance 
of finding a typical trajectory in a particular area of phase 
space, provided that it has not escaped until time t. 



B. Random maps 

Consider choosing a random sequence of maps f t by vary- 
ing the parameter a — a t — a* + 6£t- ft(xt) = f{xt\5£,t)- At 
each time t applying the same perturbation £ t to all N trajec- 
tories, £t,n = 6 m Eq. ([T]), corresponds to the random map ap- 
proach. The sequence of random perturbations can be indexed 
by the realization r as £ t ,r, with which we have different re- 
alizations of the sequence of random maps: f t r . For a fixed 
realization, we can again consider the set of points, initialized 
at a particular time t = t*, that never escape r for t ±oo. 
This set is called a snapshot saddle |[T5l . whose geometry gen- 
erally changes with time t. Its unstable manifold at time t, 
from which points never escape r for t — > — oo, can be seen 
as the open map counterpart of snapshot attractors of closed 
random maps |5 1. This way the snapshot attractor is said to be 
defined in a pullback sense as the set which is approximated 
by identically perturbed trajectories initialized in the infinitely 
distant past, and more recently it has been referred to as a pull- 
back attractor ifTTl . The constituent trajectories are distributed 
according to a time-dependent sample measure I12T1 . which is 
a generalization of the Sinai-Ruelle-Bowen (SRB) measure of 
the autonomous case (2]|2T]. The unstable manifold of the 
snapshot saddle, too, depends on the entire history of £ t r over 
t E (— oo, i*] - whereas the snapshot saddle itself depends 
also on the future t G [t«, oo). 
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For an ensemble of trajectories initialized at t = 0, later 
at some t ^> the normalized density of surviving trajec- 
tories p(x, t) will be concentrated around the unstable mani- 
fold belonging to the time t snapshot saddle, distributed ap- 
proximately according to a time- and realization-dependent 
c-density /5 t r (x). The associated time-dependent generaliza- 
tion of c-measures supported by the unstable manifold obey 
the following relation: 





autonomous maps 


random maps 


noisy maps 


attractor 


strange 


snapshot 


noisy 


measure 


SRB 

p att 


t-dep't sample- 

~att 


fuzzy 

~att / » ati\ / ,-.att\ 

V — \Mt,r/r — \Mi,r/t 


saddle 


fractal/chaotic 


snapshot 


noisy 


measure 


c-measure 

A* 


f-dep't c- 


fuzzy c- 

A = {pt,r)r = (£h,r)i 



p, t +l, r (A) 



TABLE I: Summary of measures of chaotic systems discussed in the 
(3) text. The measure in association with the saddle is supported by the 
unstable manifold of it 1151 . 



We say that the measure fi t r is conditionally-map-invariant, 
but it is not time-invariant. 



III. ESCAPE RATE 



C. Noisy maps 

The physical picture for noisy maps is provided by molecu- 
lar diffusion, in which case random perturbations act indepen- 
dently on each particle. In terms of the dynamics described 
by Eq. (JTJl, this means that t; t ,n and £, t ,n' are independent for 
any pair of n ^ n'. From the point of view of the random 
maps, the noisy map corresponds to combining the N tra- 
jectories of all R realization together, with i?, AT — >• oo. In 
case of attractors, this corresponds to combining the snap- 
shot attractors to build up the so-called fuzzy attractor lfl2l . 
and the natural measure supported by the fuzzy attractor is 
the average of those supported by the snapshot attractors: 
p att (A) = {$K{A)) r = (A?* r *(A))i El ED. This naturally 
extends to the case of open maps, where the c-measure p, of 
any set A C r of the noisy maps is given by: 

fi{A) = (JhA A ))t = (faA A ))t, (4) 

where the last equality is guaranteed by the ergodicity of £ 
and shows that p, is naturally time-invariant. This means that 
the normalized density of surviving trajectories in the noisy 
map p(x, t) converges p(x) for t — > oo, where p = (pt, r )r = 
(pt,r)t * s independent of time or realization. Note that we 
use the following notation for averaging with respect to, e.g., 
realizations: 

R 

(•)r=lim (•}*!= lim ]T. 

r=l 

A summary of the relevant measures mentioned in this sec- 
tion is given in Tab. [I] In the remainder of this paper we dis- 
cuss two fundamental quantities of the dynamics, the escape 
rate k and the dimensions D of the relevant fractal sets. We 
are mainly interested in comparing results observed in random 
maps to the corresponding noisy maps (for a fixed distribu- 
tion of £ and fixed S), and also compare these two cases to 
the unperturbed map for increasing values of the perturbation 
strength S. 



A. General relations 

A signature of chaoticity in autonomous open systems is 
the exponential decay of the survival probability inside r for 

t — >• oo: 

P(t)= lim N(t)/N(0)~exp(-Kt), (5) 

JV(0)->oo 

where n is the escape rate. We call these systems fully chaotic. 
In terms of the analysis based on surviving trajectories pro- 
posed in Sec.[n] the exponential decay in Eq. |5]) corresponds 
to a fixed fraction exp(— k) of surviving trajectories not es- 
caping after each time step. Considering that the denomina- 
tor in the right hand side of Eq. |2| is a normalization factor 
accounting for the escape of trajectories in one iteration of 
/, one obtains the well-established relation for autonomous 
maps lfT3l : 

K = -ln/i(/- 1 (r)). (6) 

In random open maps, the results of Sec. [II] show that 
Eq. (|6]l can be applied for each realization leading to a 
time- and realization-dependent single-step escape rate kt,r = 
— In fi t r (f^ (r)). In the spirit of Eq. (j^, the physically rel- 
evant escape rate for a fixed realization r is obtained aggre- 
gating the escapes over different times. Therefore, the overall 
escape rate of the random map k is given by: 

K = (k t ,r)t = («i,r)r = faAftp ( r )))r, O) 

where we used the ergodicity of the random perturbation £ as 
in Eq. (j4j). Ergodicity guarantees that all (typical) realizations 
of £ lead to the same escape rate k over t = 0, . . . , T — > oo, 
and that this value equals the mean obtained over different 
realizations r = 1, . . . , R oo. 

For noisy maps the c-measure is time-invariant and the ex- 
pression corresponding to Eq. |6]l is as follows: 

k = -\n{p t Af^r{r)))r, (8) 

where we used Eq. 
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Comparing Eqs. |7]) and ([8]l we see that the difference 
stands in the order of taking the average and the logarithm. 
We now rewrite Eq. (|7]i as the logarithm of a geometric mean: 

K = \im^ - In \J[ h,r 1 (O) J . 

which is known to be always smaller or equal than the arith- 
metic mean used for k in Eq. ([8]). With this, we arrive at our 
first result: 

k < k, (9) 

i.e., the escape rate k in the noisy-map configuration (pertur- 
bation applied independently to each trajectory) is smaller or 
equal to the escape rate k in the random-map configuration 
(perturbation applied consistently to all trajectories), for the 
same random process £ t and perturbation strength 8. In fact, 
this inequality is due to the concavity of the In function - the 
same way as the inequality of the arithmetic and geometric 
means. Equality is achieved only when fl t r (f^(r)) is in- 
dependent of time. Typically, however, for increasing 5 we 
expect k to become increasingly smaller than k. 

B. Finite-size estimation 

Important properties of finite-size S estimators include 
the bias or accuracy and the variance or precision, which are 

( S) 

respectively given by the expected value (e r ) r (minus the 

true value) and variance var[e r ] r of a distribution created 
by an ensemble of realizations of a relevant quantity. Each 
realization is produced by assigning random values to mem- 
bers of a finite-size set of the relevant quantity (e.g. initial 

(s) 

conditions, sequence of perturbations, etc.). If e r converges 
to the true value for any r as S — > oo, then the estimator is 

said to be consistent. If (er ) r equals the true value for any 
S, then the estimator is said to be unbiased. When the esti- 
mation would involve a finite R number of realizations, e.g. 
by simply taking the mean over different finite-size estimates, 
then to work out the improved precision of this estimation we 
have to consider the combined estimator e^ s,R ^ = (er )£L 1 , 
and the distribution created by an ensemble of makeups of R 
realizations each. With the standard terminology a makeup is 
then a realization of a group of realizations. 

In this subsection we consider the nontrivial case of the 
random maps only. In the previous subsection we saw that 
for every (typical) realization r the escape rate converges 
to the same value k = (kt, r )t m the limit of observation 
time T — > oo. In practice, T is restricted to a maximum value 
T max due to the finite number of initial conditions N(0), 
which, according to Eq. Q, is: 

Tmax ~ K^ 1 IniV(O). (10) 

Due to the ergodicity of £, the same k is obtained averag- 
ing the time and realization-dependent single-step escape rate 
over different realizations k = (k t: r)r- We thus see that there 



are two possible strategies to improve the precision of esti- 
mating k: (i) increasing N(0) (which, for simplicity we de- 
note hereafter by N) or (ii) increasing R. In this section we 
discuss in detail the finite-T, -R and -N estimation of k, as 
well as the precision of estimation and its scaling with T, R, 
and N. 



1. Accuracy of estimation 

It is useful to distinguish between two steps in the estima- 
tion of k: the first corresponds to the estimation of k t r , the 
escape rate for a single iterate and for a single realization of 
the random map, with a finite number of initial conditions N, 
and the second corresponds to the averaging of k t>r over a 
time interval of length T and in turn a number of R different 
realizations. 

The first step applies to the case of autonomous maps as 
well, and for simplicity we discuss this step in the frame- 
work of autonomous maps. The escape rate k is estimated 
through estimating the measure inside the escape region, ac- 
cording to Eq. After sufficiently long times, the c-measure 
i s estimated simply as the fraction of the surviving 
trajectories with one iteration of the map, i.e., N(t + l)/N(t). 
With different realizations of the (finite number) of initial con- 
ditions we expect N(t + 1) to feature a binomial distribution 
B(N(t),p(f- 1 {r))), whose mean is A/Yi)/^/- 1 ^)). This 
shows that N(t+ 1) /N(i) is an unbiased finite-iV estimator of 
/i(/ _1 (r)). However - lnN(t+l)/N(t) is a biased (inaccu- 
rate) finite-iV estimator of k because the average is performed 
after taking the logarithm - applying once more the same rea- 
soning leading to Eq. In practice, it is important to guar- 
antee that N(t + 1) is sufficiently large so that this bias is suf- 
ficiently small. Numerically this is not always easy because 
N oc exp(— Kt) and time t must also be sufficiently large in 
order for N(t + 1) /N(t) to be a good estimate of plf^ 1 (J 1 )) 
(convergence of the initial density to the c-measure p). 

Assuming we have accurate estimates of K t r , we proceed 
to the second step and consider the effect of averaging over 
time (kt, r )t or realizations (k t ^ r )r- As argued by Eq. (|7j, the 
ergodicity of £ guarantees that both averages converge to the 
same value k. In practice it is also interesting to consider finite 
averages performed simultaneously over T time steps and R 
realizations, resulting in the combined estimator: 

= <tor>K-l>£l- (ID 

This also converges to k for either T —> oo or R — > oo, which 
makes it a consistent estimator. Moreover, the finite-T and 
finite-i? estimations are also unbiased. For this we consider 
the values of k^ T,R ' obtained for a set of m — 1, . . . , M — > oo 
makeups, each one with (fixed) R realizations and T time 
steps for the estimation. According to the central limit theo- 
rem, the distribution of the estimates with respect to the make- 
ups m will be approximately normal with average equal to k. 
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2. Precision of estimation and its scaling with T, R, and N 

The variance of the same distribution scales with the num- 
ber of terms being averaged as follows: 



'(«) 



EE <(/^*)) 2 } m - (, 



TR' 



(12) 



where a\ = <t?.v (T = 1, R = 1) = var[k t , r ]t — var[k t , r ]r- 
We now compare the the two different strategies of improv- 
ing the precision of estimating k: (i) R fixed and T —} oo; and 
(ii) T fixed and R oo. For a fixed N number of trajecto- 
ries all trajectories escape in some finite time T max given by 



Eq. ( 10 1. Substituting the latter into Eq. ( 12 1 we arrive at the 



scaling law: 



o? &> oc 1/lniV. 



(13) 



The latter indicates that with strategy (i) a steady improvement 
of the precision of estimates can be achieved by increasing 
the number of trajectories N with an exponential rate. Thus, 
the precision can be improved much more effectively using 
strategy (ii), i.e., by increasing the number of realizations R, 



so that the improvement, according to Eq. ( 12 1, is (inversely) 



proportional (as opposed to a logarithmic relation) to R. 



C. Examples 

To illustrate aspects of transient chaos in fully chaotic open 
systems an area-preserving baker map will be analyzed. A 
general property of the baker map, (x t+ i, yt+i) = B(xt, yt), 
is that the mapping rule B is defined in a piece-wise manner, 
such as: 

B_ = (x t /a,ay t ), y t < 1/2, (14a) 
B+ = (1 + (a: ( - l)/o, 1 + a(y t - 1)), y t > l/2(14b) 

where a is the only free parameter. Trajectories mapped out- 
side the unit square (x, y) € [0, 1] x [0, 1] are considered to 
have escaped (open boundaries). 



Example 1: Area-preserving naturally open baker map 

In this example we obtain an analytic expressions for the 
escape rate that illustrates inequality (|9j. For a > 2 the map is 
said to be naturally open. The stochastic perturbation is added 
in a around a fixed value a* as at = a* + S£t, where a t > 2, 
and £ t is an independent identically distributed (iid) random 
variable with zero mean and finite variance <7£. 

First, consider the random map. Starting with the unit 
square, after T iterations a number of 2 T strips of equal 
width 1 ja t remain. The escape rate can be thus obtained 
as EE]: 



In — ) = ( In 
v a t / t \ 2 



--*(£'>■ (l5 ' 



where = ln(2/a*) is the unperturbed escape rate, ergodic- 
ity Q has been used, and the approximation is obtained as a 
second order Taylor expansion in 8. It is worth noting that the 
escape rate decreases with the perturbation intensity, i.e., the 
trapping is enhanced by the perturbation. A comparison of Eq. 
( fl"5j ) with Eq. |7]l reveals that the fraction of surviving trajecto- 
ries in terms of the c-measure is /i t (/ t _1 (I^)) = 2/a t . Notice 
that in this simple example the c-measure depends only on 
the current value of the perturbation, but not on the complete 
history. It is a consequence of the fact that at each time step a 
fraction 1 — 2/a t of the surviving trajectories escape. The log- 
arithm of the survival probability In P(t) can be thought of as 
a simple random walk (with a drift k) and therefore the scaling 
laws discussed in Sec. |HIB~2| can be obtained explicitly. The 
reduction of k with 8 (noise-enhanced trapping) can be under- 
stood in this case simply as a consequence of the concavity of 
the logarithmic function n = ln(a/2). 

Next we consider the noisy map, apply Eq. (|8j, and take the 
following approximations: 



= - In 



-In 



<5(7£ 

a* 




(16) 



Interestingly, for small perturbations the noise increases the 
trapping by reducing k* by twice the amount as in case of 
random maps. Choosing £ to be uniformly distributed in 
[—1,1], we can also compute k and k exactly using (. . .) t = 
(. . .),. = g J j . . . dl;. For a* = 2.5 and S = 0.5 we obtain 
k = 0.2164... and k = 0.2096.... 



Example 2: Area-preserving closed baker map with a leak 

In our second example we explore a case in which no sim- 
ple analytic expressions for fj, tj r and k exist, but in which nu- 
merical results confirm the validity of the scaling <j 2 k) ozl/T 



(12 1, the inequality k < k d9j, and noise-enhanced trapping 
for small <5. We start from the area-preserving closed baker 



map, obtained by fixing a = 2 in Eq. ( 14 1, into which we 



introduce a leak of width A x =0.1 vertically in the mid- 
dle of the map's phase space. Any trajectory that is mapped 
into the leak region is considered to have escaped. Differently 
from the previous example, here the stochastic perturbation 
is defined to act on the coordinates independently, such that: 
x ->• x + 8£ x and y -» y + 6£ v , where £ x ,£ y e [-1,1] 
are uniformly distributed independent random variables. This 
way escape out of the unit square is also possible. 

In Fig. [T] (a) we show the temporal decay of the number 
of surviving trajectories N(t) for R — 10 4 different realiza- 
tions of £ t (all iV(0) = 10 6 trajectories are exposed to the 
same sequence £t) ll24ll . In agreement with Eq. the es- 
cape rate of the random map k can be obtained by averag- 
ing the logarithm of the number of surviving trajectories over 
different realizations, i.e., (In N(t)) r ~ lniV(0) — kt. In 
Fig.[T](a) this slope is shown to be somewhat larger than the 
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value R which in principle (for N(0) — »• oo) corresponds to 
independent perturbations applied to each initial condition, so 
that \n(N(t)) r ~ ln7V(0) - kt. It is expected that k is ob- 
tained here rather inaccurately, because our procedure acts as 
if k was determined from a small R = 10 4 number of trajec- 
tories, and the value was redundantly replicated-and-averaged 
iV(0) = 10 6 times. Nevertheless, the results are suitable to in- 
dicate the general inequality ^9), In the insets of Fig. [T](a) we 
explore the dispersion of In N(t), which is shown to follow a 
t 1 / 2 scaling, according to the general rule ( 12 1 and observing 
Eq. Q. Notice that proper scaling is attained after a rela- 



1 1 1), once the initially 



tively short transient time (£' as in Eq. 
uniform densities closely approximate the c-density. 

The dependence of k (□) and k (•) on the noise strength 
6 is shown in Fig. [T] (b). The difference between k and 
k steadily increases with the perturbation strength and both 
curves show a nonmonotonic dependence with a local min- 
imum. This effect has been explained in [20 1 as follows: 
for small 5 the probability of escape through the leak is re- 
duced k > Kg=o because the invariant density over the leak is 
smoothed in comparison with the unperturbed case which has 
a fractal support; for large noise more and more trajectories 
escape through the open boundaries, leading to an increase in 
k. Here we show that k follows the same general trends, but 
the noise-enhanced trapping is somewhat less effective due to 
the inequality (|9j. In the next section we explore the effect of 
the random perturbation on the associated fractal dimension. 



IV. DIMENSIONS 

A. General relations 

An important difference between the dimension and escape 
rate is that the noisy perturbations wash out the fine details 
of the invariant sets and the (asymptotic) dimensions are not 



fractional D 



D. 



phase space 



Speaking about fractality 



is thus meaningful only in terms of the random map, and a 
nontrivial inequality like Q cannot be established for the di- 
mensions. Here we discuss three different relationships be- 
tween the dimension and other quantities. In the next section, 
subsequently, estimations of the dimensions will be based on 
these relationships: (I) based on box-counting (BC), (II) the 
Kantz-Grassberger (KG) relationship, and (III) based on the 
uncertainty exponent (a). 

(I:BC) The fractal dimension is defined directly by the rele- 
vant measures discussed in Sec. |ll] The full fractal dimension 
spectrum D q of a measure can be computed by applying e.g. 
a box-counting algorithm. The relevant measures in case of 
random open maps are the sample measure supported by the 
random saddle and the c-measure supported by the unstable 
manifold of the random saddle. The dimension spectrum of 
the sample measure in the stable direction is identical with 

(2) 

that of the c-measure, denoted as D q ' . Although large scale 
features of fit in phase space, in association with a limited 
memory of the past, depend on time, its asymptotically fine 
details are determined by the complete history of the system's 
evolution. The fractal dimension quantifying the asymptotic 



scaling of the measure is thus constant fill . This too makes 
the dimension a qualitatively different characteristic number 
from that of the escape rate k t , which changes in time together 
with/t t [Eq. Q]. 

(II: KG) For autonomous fully chaotic 2D open maps the 
connection between the geometry and the dynamics on the 
chaotic saddle is established by the pair of the Kantz- 
Grassberger relations, which reads as follows [26 |: 



Df» = 1 - 



A' 
A — K 



(17a) 
(17b) 



where D^p and are the partial information dimensions 
across (along) the stable (unstable) and unstable (stable) man- 
ifolds, respectively, and A and A' are the corresponding aver- 
age positive and negative Lyapunov exponents on the chaotic 
saddle, respectively f27l . In conservative area-preserving sys- 



tems A = —A', as a result of which: — D[ z> , and e.g 



l(2) 



Dq 1 ' = D^' for any q. In case of random maps, similar rela- 



tions to ([T7]) hold f6): 



D 



(i) 



D 



(2) 



A 



K 



|A'| 



(18a) 
(18b) 



Analogous to kt, r (except for 12711 ). we can define the time- 
and realization-dependent one-step positive and negative av- 
erage Lyapunov exponents as ensemble averages, e.g. A t , r — 
(\t, r ,n)n (n as in Eq. ([IJ), and with this A = (Xt. r )t = 
(^t,r)r, where the latter equality is due to the ergodicity of 
noise. The same applies to A'. For each trajectory At r n = 

lny t , r , n and \' t r n = ln$ rn are defined by y tyr ,n and y' t>r<n , 
which are respectively the stepwise stretching and shrinking 
rates along the corresponding covariant Lyapunov vectors, 
i.e., the corresponding manifolds. 

(Ill: a) Fractality can be related with the concept of uncer- 
tainty. The latter is measured by the uncertainty exponent a, 
which specifies the scaling of the ratio of the number of un- 
certain boxes to that of all the boxes with the resolution or box 
size, such as: Nb{s) / N$(e) oc e a . The certainty of a box is 
defined so that any trajectory from it takes the same route of 
escape. For example in case of the unperturbed baker map, 
Example 1 of Sec. IIIC trajectories can escape from either 



the left or the right side of the leak. The uncertain boxes for 
e — > shrink onto the stable manifold of the chaotic saddle, 
and thus the uncertainty exponent is related to the Hausdorff 
dimension as |[T6l : 



D 



1 



(19) 



because the number of all boxes in a plane scales as: Nq{e) oc 
e~ 2 . The uncertainty grows with decreasing a, which means 
that it gets more difficult to improve the predictability of the 
outcome for the typical (randomly chosen) trajectory by in- 
creasing the precision of the initial condition. Equation ( 19 1 
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FIG. 1: (Color online) Difference between the escape rate in the noisy (k) and random maps (k). (a) Number N(t) of surviving trajectories 
out of N(0) — 10 6 over time for perturbation strength 5 = 0.1 in the area-preserving closed random baker map with a central vertical leak 
of width A x = 0.1. A number of R — 10 4 experiments (realizations of f) were carried out and (200 of them) are shown as gray lines in the 
backdrop. The lines in the front show: (jV}^ =1 (blue solid line), corresponding to the noisy map as in Eq. 181; exp (In N}^ =1 (red dashed 
line), corresponding to the random map as in Eq. I7J; and exp((ln N}^ =1 ± a\ n jv/2) (dash-dot lines), corresponding to the standard deviation 
around the expected value in the random map case. Lower inset: histogram of N(t — 15); upper inset: scaling of the dispersion criniv (half 
the vertical space between the dash-dot lines) as t 1 ^' 2 (the solid straight line indicates a slope of 1/2). (b) Dependence of k and k on the 
perturbation strength S. A number of R — 10 4 x 6 experiments are done for each 8 > (and this time for the noisy map all N(0) = 10 6 
trajectories are perturbed independently, in all the R experiments). The escape rates k and k were estimated from data from time t = 6 up to 
20 (giving T — 14). In case of the random map errorbars indicate the standard deviation ov^n . In case of the noisy map the errorbars a^) are 
smaller than the square marker. In (a) and (b) data points are connected by lines to guide the eye. 



indicates that the stronger the fractality the greater the uncer- 
tainty, and the same relation holds for noisy as well as random 
maps. For noisy maps predictability cannot be arbitrarily im- 
proved by improving the precision in the initial conditions, 



consistent with a — and D, 



n(i) 



1, If the perturbation his- 



tory is known, then the random map framework is relevant (in 
which case the term 'random' is rather misleading), fractality 
is resolved (Dq < 1), and predictability can be improved 
(a > 0), similarly as in the unperturbed case. 



B. Finite-size estimations 



Different algorithms for the computation or estimation of 
the fractal dimension based on the relations of the previous 
section (I,II,III) are commonly used. Here we discuss the 
accuracy (bias and consistency) and precision (variance or 
spread) of the estimations based on these algorithms when ap- 
plied to random maps. As discussed before [see Eq. (10l], 



with finite N number of trajectories numerical experiments 
are limited to a maximum time T max ESIl . and estimations 
are based preferably on averaging over different experiments, 
that is, realizations R of the perturbations, labeled as strategy 
(ii) in Sec. [mBl 



1. List of estimators 

(I:BC) Most naturally the fractal dimension is calculated by 
a direct estimation of the scaling of lxiNb(e) with — me (or 
that of the information for the information dimension). A/j, is 
approximated by the count of boxes in a regular rectangular 
grid that contain at least one point out of finite N(T < T max ) 
points that represent e.g. the unstable manifold. The scal- 
ing line is fitted by a straight line over a finite e-range, up to 
a minimal box size whose slope estimates the dimension 

D^Cr ■ acma l estimator is then defined as the average 
over R realizations: 



Dbc = D 



BC 



BC,r /r=l 



(20) 



(ILKG.a) Another type of estimator is based on the KG re- 
lations ( |18) , For simplicity, we focus only on the second of 
these relations (the same conclusions apply to the other op- 
tion). Similarly as for Dbc, we define the estimator as an 
average of the single-realization KG estimator over R realiza- 
tions: 

/ C (T) s (T) \ R 
A = A(2,T,_R) _ / A(2,T) ,R _ / A r - K r \ 

^KGa — ^KGa — \ U KGa,rl r=l — \ |j\'( T )| / ' 

r (2i) 

where k[ T ' = (kt, r } t t J~J' (t' > 1 for well-approximating 
the measure) is a version of k^ T ' R ' 1 of Eq. ( jll] ), without 
averaging with respect to the realization, and similarly e.g. 



{(T) 
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(ILKG.b) An alternative KG estimator could be defined by 
changing the order of division and averaging: 



D 



KGb 



D 



(2,T,R) 
KGb 



(A 



(T) 



Kr ) r —\ 



\\>(.T,R)\ 



<ia; (t) i>? =1 



(22) 

where the last equality follows from the convention ( [TT] ). Note 
that here we do not have a single-realization estimator for the 
dimension, but only for the escape rate and the Lyapunov ex- 
ponents. 

(Ill: a) The stable manifold is contained by the basin bound- 
ary, which separates regions of initial conditions from which 
trajectories take different routes of escape. Thus, the num- 
ber of uncertain boxes can be approximated also by an ap- 
proximate finite-e m j„ resolution survey of the basin bound- 
ary. Trajectories are initialized on a regular rectangular array 
of points, and all of them will escape one way or another in a 
finite time. The route of escape will assign logical values to 
the grid points. The scale is varied by a rectangular grouping 
of k 2 grid points, so that e = ke m % n . A box of size e is certain 
if all the k 2 grid points are assigned the same logical value. 
This way we have another means of box-counting estimation 
of either the uncertainty exponent, or, according to relation 
dl9), the Hausdorff dimension of the stable manifold: 



)r=l- 



(23) 



2. Accuracy of estimation 

Similarly to the discussion of the escape rate, we first 
discuss the accuracy of the finite-time estimation, and then 
the accuracy in relation with averaging over realizations. A 
single-step and single-realization measure of the dimension, 
akin to k tr , has not been defined, because the true dimension 
is constant. However, the single-realization box-counting es- 
timate Dg G *J is found to be time- and realization-dependent. 
This can be explained with its relationship to the finite-time 

" (2 T) 

single-realization estimator D KGci as follows. First, in terms 
of the baker map, which is a paradigmatic example for e.g. the 
fractal structure of chaotic attractors or saddles, a number of 
T iterations determine the scaling of ge ometry down to a size 
e oc |A'| T (see Example 2 of Sec. IVC I. Second, with a finite 



N number of points to evaluate D^F' , fractal scaling holds 
down to a scale related to N as: In 1 je„ oc In N. These two 
points suggest that with finite N, the system has a memory of 



T oc lnl/e, oc IniV 



(24) 



steps back into the past. Indeed, we have found numerically 

~ (2 e ) A (2 T) 

^BCr t0 correlate with D y KG( ^ r for an appropriate T. Thus, 
the time-dependence of the latter can be related to the time- 
dependent finite-memory/scale properties of the measure. In 
the limit N —> oo, the entire past determines the dimension, 

" (2 e ) " (2 T) (2) 

and both D BC *' and D KG ^ converge to the true value D\ 
as Inl/e* and T — > oo, for any realization. For finite T and 



- (2 T) it) 

N, the accuracy of D KG ^ depends on the accuracy of k r 



" (T) 

and \r .It has been argued in Sec. 



IIIB 1 



(T) 

that kr is biased, 



but it can be effectively reduced by choosing N large enough. 
The situation with Xi^ is similar, and therefore, so it is with 

* (2 T) A (2 s ) 

^KGar- As f° r "gQr ' some box-counting algorithms has 
been reported in QUI to be negatively biased, which agrees 



(2,T) 
KGa,r> 



with our numerical experiences. Unlike the bias of D 
the bias of Dg G "J may not be insignificant relative to other 
components of the bias of the total estimator Dg G " R \ 

Fitting the approximate scaling line of A4(e)/-Vo( £ ) to ca l- 
culate the dimension through the uncertainty exponent, ac- 
cording to Eq. p3] >, leads to a different bias than what we 
have if directly the scaling line of JVb{e) is fitted. 

Next, we focus on the second step, i.e., the averaging over 
different realizations. We start with com paring D KGn and 

that k (T > R) 
N -> oo). 
Therefore, 



IIIB 1 



DuGb- First we recall our results from Sec. 
is an unbiased estimator of k (assuming that 
We can argue that the same applies to A^^. 
DjiGb — > D as R — > oo, and we conclude that L>KGb is a con- 
sistent estimator. Now we notice that Z?a'G6 an d D^Ga will 
typically lead to different values because the ratio of the aver- 
ages is different from the average of the ratios. For the case of 
iid random variables we have: (x/y) = (x)(l/y) > (x)/(y). 
(The inequality holds also for dependent x and y.) Therefore 
we conclude that Dkg<i is not a consistent estimator, and it 
always overestimates the true value. The bias of these esti- 
mators at finite R can be determined by considering the aver- 
age value of Dxca an d D^Gb over different makeups (as in 



Sec. Ill B 1 1. In this case, both estimators would correspond 
to the average of ratios instead of the ratio of averages. There- 
fore, we conclude that both Dkgo, an d DxGb a re biased esti- 
mators. 

As for Dec, even if we do not have an analytical model of 



it, the correlation between Dg G *J and D^ G ' a suggests that 

Dbc should also be expected to be biased (beside the effect 
reported in [30]). Because of its similar nature, the same can 
be said about D a . 

In the example above the source of the bias was taking 
the arithmetic mean (. . and it is natural to ask our- 

selves whether other approaches could lead to a better esti- 
mate. In this regard e.g. the harmonic mean, 
is not expected to perform better when the numerator is 
not constant. In case of iid random variables x and y, 
again with a reference to the numerator and denominator of 



}(2,T) 



D 



(2,T) 



the harmonic mean underestimates the true value: 

^ l„\ Hf|, „ J //J2 



KGa,r> 

(y/x)- 1 = (y)- 1 ^- 1 )- 1 < (x)/(y). When o*/(x)- is 
much greater (smaller) than a 2 /(y) 2 , the arithmetic (har- 
monic) mean gives the better approximation, since we find 
that e.g. (y) / (y^ 1 )^ 1 — 1 oc a 2 /(y) 2 . It is easy to verify that 
for equal values of these, the harmonic mean yields smaller 
bias. 

* (2 T R) 

To the end of creating an unbiased D y KG ' a -type estimator, 
an appropriate generalized /-mean can be found only when 
the numerator and denominator are interrelated in a special 
way. As an example, when the Jacobian \J t \ = yty' t of a 
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closed dissipative system is constant, as in case of the fluid 
flow in [5 1, the appropriate generalized /-mean can be char- 



acterized by the form: f(D^ 



D 



(2,T) 



/{D 



(2,T) 



1). In 



the latter D 



(2,T) 



" (2 T) 

tained from D K ' G ^ r when n 



is the Lyapunov dimension that can be ob- 
0. 



3. Precision of estimation and its scaling with R and N 



" (T) (T) 

Considering that e.g. X r is defined similarly to k r , its 



variance scales similarly as given by Eq. (12 1, which is in 



herited, along with the Gaussian form of limit distributions, 

~ (2 T) 

by DjQGa r- Considering averaging with respect to realization 
too, we have: 



2 J_ 



(25) 



For some fixed R, using relationship ( 24 1, we have then: 



a 2 h <x lnl/ff* oc 1/lniV. 



We note first that these scaling laws of precision apply to at- 
tractors of closed systems too, in agreement with previous 
derivations by Namenson et al. [5|. Second, they are formally 
analogous to scaling laws of the escape rate estimates, Eqs. 



( 12 1 and ( 13 i, respectively. Therefore, third, just like with the 
escape rate, it would be an overwhelming numerical burden to 
improve the precision of the dimension estimate by increasing 
the ensemble size. It is done much more effectively by pro- 
ducing a number of estimates with different realizations of the 
perturbation sequence and averaging them. However, unlike 
in case of the escape rate, this approach for estimating the 
dimension introduces a bias, as discussed previously in Sec. 
IIVB II 



C. Examples 



We consider here the two examples specified in Sec. Ill C 



Example 1: Area-preserving naturally open baker map 



prompted as follows. (For simplicity R = 1, and we omit the 
realization index r in the notation.) The width of the stripes 
can be considered as the minimal box size to cover the order- 
T approximant of the manifold: = XlJ =l l/at- Thus, not 

r (2 T) 

only the dimension estimate D ' ' depends on the realiza- 
tion of the random sequence at, but also the minimal box size 
e*. Therefore, it is not straightforward to generate a distri- 

" (2 T) 

bution of Dq belonging to the same just by generat- 
ing a large M number of 'realizations'. However, by increas- 
ing the sequence length T, the effects of decreasing box size 



can be explored. In Fig. 



| (a) 



dimension estimates D, 



(2,T) 



are plotted against corresponding minimal box sizes, or rather 
In For each fixed T the data points align to a hyperbola 
segment. The large number of realizations generate sampling 
distributions of Dq ' and or In . The vertical dis- 
tance between a pair of dash-dot lines indicate the standard 

" (2 T) 

deviation of the sampling distribution of Dq ' . The upper 
inset in the figure shows the standard deviation df )o against 
In belonging to fixed T's, and confirms that scaling law 



(26) (26 1 holds. The lower inset shows the sampling distribution 



of D^-°^>, which can be well-approximated by a Gaus- 
sian form. 

We can obtain the scaling laws formally as well, as fol- 
lows. We shall introduce now the following notation involv- 
ing the minimum box size: v = lnl/e*, and for the baker 
map we have that: v = T(\n a t ). It follows then that: 
ct 2 cx T 2 /T = T. Next, Eq. (|27]i can be rewritten as: 
Do = ln2T/iA If we linearize this about the mean D, we 
can establish that: <j 2 - 

JJo 

T limit: cr 2 - 

Do 

have that D oc T, and with this: tx \/D, which conforms 



\D' ^(D)\ 2 o-l. That is, in the large- 
\n z 2T 2 /D i T. From the definition of v we 



with Eqs. (|25j and ((26 1. 

We note that in this special case, considering Eq. ( |27| ), the 



harmonic mean of the D 



(2,T) 
KGa,r 



(or Dg£*J) values would pro- 



" (2) A (21 

vide an (approximately) unbiased estimator of D\ = D Q . 
A comparison of the harmonic and arithmetic means is shown 
in Fig. |2](b), where their difference for decreasing T or In l/s 
is clearly indicated, which gives the measure of the bias of the 
arithmetic mean. 



By this example we intend to support the claim that dimen- 
sion estimates are generally biased, and also the scaling laws 
of the precision of estimates given by Eqs. ( |24] > and ( |26] i. The 
finite-time positive Lyapunov exponent can be obtained such 



A 



(T) 



(In a t . 



/*=!• 



Note that after T iterations each 



Example 2: Area-preserving closed baker map with a leak 



The fractal dimension D 



(i) 
o 



" (2) 

Dq ' is evaluated for the 



stripe is stretched by a factor of Hf =1 at- On the other hand, 

(T) 

having already obtained an analytic expression for given fined in Sec [TvbT] distinguished by different markers in Fig 



same values of the perturbation strength S as considered pre- 
viously in case of evaluating the escape rate (Fig. (T](b)). It has 
been estimated in five different ways using the estimators de- 



by Eq. ( 15 1, we can write the order-T approximant to the di- 



mension across (along) the unstable (stable) manifold as: 



A(2,T) _ A(2,T) 



In 2 



(27) 



We numerically generate a large M number of sequences of a t 



of length T, by which the first part of scaling law ( 26 1 can be 



[3] (The reader is referred to the Appendix for details.) The re- 
sults indicate a characteristic enhancement of fractality, with 
a maximal value of D for some finite 5. This value should be 
the same as the one for which maximal trapping occurs, for 
the following reason. With additive perturbation of the baker 
map, the Lyapunov exponents are unchanged, which is not 
simply because the perturbation terms do not explicitly appear 
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200 300 



FIG. 2: Fractal dimension in the area-preserving naturally open baker 
map. (a) Fractal dimension Z)q 2,t ' vs minimal box size e*. A 
number of M — 10 4 experiments (realizations of a t ) were car- 
ried out to generate statistics. The standard deviation of _Dq 2 ' t ' 
with T fixed is marked by a pair of dashed lines. Lower inset: his- 
togram of Z)q 2 ' t ~ 500 '; upper inset: scaling of the dispersion ajj 
(half the vertical space between the dash-dot lines, approximately) 
as (In l/e*) _1//2 for large lnl/e*'s, the e*'s belonging to the re- 
spective expected values [if )g (the solid line has a slope of 1/2). (b) 
Comparison of the harmonic (lower line) and arithmetic mean (upper 
line) estimates. For both panels T = 10/, 2 = 1,..., 50. Discrete 
data points which indicate the means and standard deviations of di- 
mension estimates are connected to guide the eye. 



in the Jacobian matrix, but also because the Jacobian matrix 
is the same in every point of phase space, for which reason 
the randomly perturbed trajectories do not have an influence 
either. Thus, considering the KG relations (fTSl) for random 



Thus, considering the KG relations 
maps, the dependence of the dimension on the noise strength 
is inherited solely from that of the escape rate in this exam- 
ple. If the Lyapunov exponents monotonically depend on the 
noise intensity, which is believed to be the case in general, 
then noise-enhanced fractality and uncertainty should always 
accompany noise-enhanced trapping. 

The different estimators seem to be biased to different de- 
grees, as the errorbars do not overlap systematically. The esti- 

A (2 T R) 

mator DjcqI (marker x ) has been shown to be unbiased and 
thus it is expected to yield the most accurate figures, closely 
approaching the true value. Therefore, its mismatch with the 
other markers indicates the bias of the corresponding estima- 
tors. Since they group fairly closely around the true value, 
they are also capable of robustly indicating the nonmonotonic 
behavior. 




0.02 0.04 0.06 0.08 0.1 



FIG. 3: (Color online) Fractal dimensions in the area-preserving 
closed baker map with a leak for several values of the additive per- 
turbation strength 5. Errorbars indicate the standard deviation of the 
sample mean in the spirit of Eq. (25), An increasing number of ex- 



periments were done for increasing 5 (R = 5 • 10 x 5, and R = 10 
for 5 = 0). Five different estimators are employed as described in the 
Appendix. All estimates are obtained for six values of 5 £ [0, 0.1] 
with 0.02 increments, but results using the different estimators for the 
same 5 are plotted with a spacing for better visibility. A gray curve 
corresponding to the best fitting quadratic polynomial emphasizes 
the nonmonotonic dependence of the dimension on the perturbation 
strength. 



V. CONCLUSIONS 

We have investigated how the characteristic measures of 
transiently chaotic systems (e.g., k, Do, A) depend on the type 
(noisy or random maps) and strength of the stochastic per- 
turbations. Random maps are described in terms of a time- 
dependent conditionally-map-invariant measure, a measure 
which we introduced as a natural combination of concepts 
from transient chaos theory lfl4l [L5l and (dissipative) ran- 
dom maps fl2l[T7ll . For any fixed time, this measure exhibits 
a clear fractal character, similarly as the measure of the au- 
tonomous system. As in the case of attractors lfT2ll . we argue 
that the measure of the noisy map corresponds to the aver- 
age of the time-dependent random map measure over differ- 
ent times (or realizations) and is smooth on fine enough scales. 
Based on this description we showed that the escape rate for 
the random map k is always larger than the one of the noisy 
map k (for the same perturbation strength). 

All measurements in numerical and experimental situations 
are limited to finite numbers of trajectories N and realiza- 
tions R. We have shown that in the random map the precision 
of finite-time estimates of k and D alike converge extremely 
slowly with N, typically as oc 1/ In N, We have shown that 
this limitation can be compensated by estimating through av- 
eraging over different realizations R, whereby the precision 
scales as oc 1/R. In the case of the fractal dimension D, how- 
ever, our results indicate that the different finite-size estima- 
tors are typically biased, and even inconsistent as R — s- oo. 
This means that the quality of estimating D, depending on 
N and R, has to be carefully analyzed, to guarantee that the 
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inaccuracy (bias) and precision (spread) are small. 

Our results regarding the dimension gain a physical mean- 
ing through its relation with uncertainty: the greater the di- 
mension, the greater (smaller) the uncertainty (exponent a), as 



can be seen by the relation a = 1 — D in Eq. ( 19 1. The finite- 
size/resolution estimates of D vary with the realization of the 
stochastic perturbation process, as a consequence of which the 
uncertainty (due to the uncertain choice of the initial condi- 
tion) depends on the specifics of the current perturbation too. 

Similarly to noise-enhanced trapping EDI , in random maps 
we have observed a nonmonotonic dependence of k and D 
on the perturbation strength 5, when the extrema occur at ap- 
proximately the same finite value of S. This entails that for the 
same perturbation strength the uncertainty is also maximal. 
The intuitive picture for this is that perturbation-enhanced 
trapping increases the chaotic life time of the trajectory, which 
makes it less predictable. The perturbation strength have to be 
increased beyond this point to steadily improve predictability. 
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Appendix 



Numerical computation of the D estimators 



Results distinguished by different markers in Fig. [3] cor- 



respond to the different algorithms described in Sec. |IV B 1 
and further details are provided below. Before each descrip- 
tion, the code name of the algorithm is followed by the type 
of marker used. 

(I:BC.a: A and y) ^ n m i s special case of the area- 
preserving closed baker map with a vertical leak centered 



around x = 1/2, the stable and unstable manifolds in 2D 
are aligned to straight horizontal and vertical lines, respec- 
tively. This allows for the collapse of all data points and box- 
counting in ID. The patterns of both the stable and unstable 
manifolds are irregular when perturbations are present, mean- 
ing that different lines that constitute the manifold may have 
different lengths. With collapsing the unstable manifold, de- 
spite its irregularity as just described, the true value of the 
fractal dimension is not affected. Numerical results, however, 
reject equality of the dimension estimates for two values of the 
perturbation strength (S = 0, 0.02). This may be because the 
irregular geometry of the manifolds introduces an additional 
bias of the very simple box-counting dimension estimator. For 
the simulation a number of N = 5 • 10 6 trajectories are ini- 
tially uniformly distributed in (x , yo) € [0, 1] x [0, 1]. The 
unstable manifolds are approximated by the trajectories which 
did not escape after tt — 20 iterations. 

(LBC.b: *) Fractality across the stable manifold can be re- 
solved also by initializing the ensemble (N = 5 • 10 5 ) along a 
line across the manifold (e.g. xq = 0.3 and yo £ [0, 1]). 

(ILKG.b: x) For the estimator D K K W ' the escape rate 



is evaluated based on the interval t G [6, tf/2], where 
t f is the minimum of all simulation run times at which all tra- 
jectories are already escaped (same simulation as for HI: a). 

(Ill: a. a: •) The basin boundary is determined along the 
same line as for (LBC.b), and the dimension of it is estimated 
by fitting the scaling line Afb(s). For this the finest resolution 
is facilitated by N = 5 • 10 5 trajectories. In case of the un- 
perturbed baker map the same basin boundary is obtained by 
checking whether the left or right side of the leak-, or whether 
the lower (y < b) or upper (y > b) regime the trajectories es- 
cape from. For the results presented the first option was taken. 
This is the estimator that overall best conforms with the sup- 
posedly most accurate estimation (ILKG.b). 

(IILa.b: +) The dimension can be obtained also by first 
estimating the uncertainty exponent by fitting the scaling line 
A/b(e)/Ao(e) and then applying relation (19 1. Interestingly, 



this approach modifies both the accuracy (leading to overesti- 
mation) and precision of estimation. 
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